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Consider a particle which is randomly accelerated by Gaussian white noise on the 
line segment < x < 1 and is absorbed as soon as it reaches x = or x = 1. 
The mean absorption time T(x,v), where x and v denote the initial position and 
velocity, was calculated exactly by Masoliver and Porra in 1995. We consider a more 
general boundary condition. On arriving at either boundary, the particle is absorbed 
with probability 1 — p and reflected with probability p. The reflections are inelastic, 
with coefficient of restitution r. With exact analytical and numerical methods and 
simulations, we study the mean absorption time as a function of p and r. 
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I. INTRODUCTION 

Consider a randomly accelerated particle which moves on the half line x > with an 
absorbing boundary at x — 0. The motion is governed by 

^=v(t), (v(t)v(t')) = 28(t-t'), (1) 

where 7](t) is uncorrected white noise with zero mean. The probability Q(x,v;t) that a 
particle with initial position x and initial velocity v has not yet been absorbed after a time 
t has been studied by McKean [l| and others DDD . For a particle initially at x = with 



Q(0,v;t) ~ ( ^- ) , t->oc. (2) 



velocity w > 0, the probability decays as 

Following Cornell, Swift, and Bray we consider a more general boundary condition. 
On arriving at the boundary the particle is absorbed with probability 1 — p and reflected 
with probability p. The velocities of the particle just after and before reflection are related 
by v f = —rvi. Here r is the coefficient of restitution, and r = 1 and r < 1 correspond to 
elastic and inelastic collisions, respectively. 

As physical motivation, we note that temporal or spatial irregularities at the boundaries 
of a system or a statistical capture process may give rise to partial absorption. The properties 
of a particle which is subject to a random force and collides inelastically are of interest in 
connection with driven granular media. 

For the partially absorbing, inelastic boundary condition the probability that the particle 
has not yet been absorbed decays as 

(2 \ 0(Pi r ) 
yj , t^oo. (3) 

Burkhardt [7| and De Smedt et al. [8] showed that the persistence exponent <p(p, r) is 
non-universal, depending on p and r according to the relation 

pr^ = 2 sin | — — — — — 7r J . (4) 
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For p = 0, 0(0, r) = i, so Eqs. (0) and © are consistent. Solving Eq. (JH) for graphically, 
one finds a unique solution < 4>{p,r) < | in the physical region 0<r<l,0<p<l, 
representing a decay slower than in Eq. (0). 
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The decay in Eq. (jSJ) is so slow that the mean absorption time 

dQ(0,v;t) 



T{Q,v) 



dtt 



dt 



dtQ(0,v;t) 



(5) 



is infinite, even for p = 0. On the other hand, the mean absorption time T*(0, v ), averaged 
over a long but finite time t* instead of an infinite time, is a well-defined quantity, which 
varies as 

V ¥ ) , t*»v 2 (6) 

and diverges in the limit if —>■ oo. 

If the particle moves on the finite line < x < 1 and is absorbed the first time it reaches 
either x = or x = 1, the probability that it has not yet been absorbed after a time t 
decays as e~ xt for long times jsj], much more rapidly than the power law (j2J) for a single 
absorbing boundary. The mean absorption time T(x,v) for a particle with initial postion 
x and initial velocity v is well defined. Masoliver and Porra 10| calculated this quantity 
exactly, by solving the inhomogeneous Fokker-Planck equation 



d d 2 



with the boundary conditions 



dx dv^ 



T(x, v) 



T(x, v) 



Til — x, —v) 



-1 



T(0,v) = 0, v<0, 



(7) 

(8) 
(9) 



corresponding to reflection symmetry and the immediate absorption of a particle with initial 
conditions x = 0, v < 0. Their result, which is rederived in the Appendix, is given by 



Ti(x, v) = Ti(l - x, -v) = i(vu; 
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(10) 



(y — x) 2 / 3 

for v > 0, where 2-^1(0, 6; c; 2) is a standard hypergeometric function 12]. The subscript 1 of 
Ti(x, v) is a reminder that the particle is absorbed the first time it reaches either boundary. 
For use below we note the asymptotic forms 

( 3 



7T 



1/2 



,1/2 







Ti(0,v) ~ ^ 



(11) 



00 
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The second of these forms is intuitively obvious, corresponding to ballistic propagation from 
x = to x = 1 in the time v~ l . 

In this paper we study the mean absorption time of a randomly-accelerated particle 
moving on the finite line < x < 1 for the more general boundary condition described 
above: absorption at the boundary with probability 1 — p, reflection with probability p 
and coefficient of restitution r. For this boundary the mean absorption time also satisfies 
the Fokker-Planck equation ((7j) with reflection symmetry (jHJ), but the absorbing boundary 
condition (jHJ) is replaced by 

T(0, -v) =pT(0,rv) , v>0. (12) 

The absorbing boundary condition © considered by Masoliver and Porra is unusual in 
that T(0,v) is specified for v < but not v > 0. That the Fokker-Planck equation with 
this boundary condition and with reflection symmetry (jHJ) is a well-posed boundary value 
problem was known long before the exact solution of Masoliver and Porra. We are 
unaware of a similar proof for the more general boundary condition (fT2|) . In this paper we 
show how to solve the more general case with an extension of the approach of Masoliver and 
Porra. 



II. INTEGRAL EQUATION FOR T(x, v) 

In the Appendix we show that the Fokker-Planck equation (J7|) with reflection symmetry 
(JH1) has the exact Green's function solution 

POO 

T(x,—v) — Ti(x,—v)+ / duuG(x,v,u)T(0, —u) , (13) 

Jo 

for v > 0, where T\ is the same function as in Eq. (FTTH) . 

G(w) = __ e ( >/ 9/ _ 1/3 (__ 

1 f X e -v 3 /9(x-y) 

~ 3 i/3r(|) I dy ( x - y y/ 3 [R{ - V) u) ~ R{1 ~ y ' u)] ' (14) 

^ _ 1 M i/2 e -" 3 /% J 15 u 3 (l-y) \ 

^' M) - 3 5/6r(i)r(|) 2/ 7/6(i_ yy/s'^y 6'6'_9y )' 1 ; 

and iFi(a,b, z) is a standard confluent hypergeometric function 
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To calculate T(x,—v) from Eq. (|13jh one must first determine the unknown function 
T(0, —u) on the right-hand side. Setting x — 1 in Eqs. (fT3 |) -(fT5 |) and substituting T(l, — 1>) = 
T(0, u) and T(0, — w) = pT(0, ru), which follow from Eqs. (jBJ) and (fT2"|) . leads to the integral 
equation 

POO 

T(0,v) ^T 1 (0,v)+p duuG(l,v,u)T(0,ru) , v>0, (16) 

Jo 



where 



G {1,V,U) = -U l/2 M l/2 e -(« 3 + U 3)/9 
07T 



9 „ / 5 7 f 3 w 3 



+ 6xF 2 1;-,-; 



(17) 



_W 3 +W 3 _ V 6 6 81 

and iF 2 (a; b, c; z) is a generalized hypergeometric function 12]. Our analytical and numerical 
predictions for the mean absorption time are derived directly from Eqs. (fTT)|) and (EI). The 
same symmetric kernel G(l,v,u) = G(l,u,v) plays a central role in Refs. |l<l llJ|. in 
determining the equilibrium distribution function P(x, v) of a randomly accelerated particle 
moving between inelastic walls at x = and x — 1 with which it collides inelastically. 



As discussed in |l4j, the quantity G(l,v, u) generalizes McKean's result for the velocity 
distribution at first return [15j from the half-line to the finite line. The probability that a 
randomly accelerated particle, which leaves x = with velocity v > 0, arrives with speed 
between u and u + du the next time it reaches either x = or x = 1 is given by uG(l, v, u)du, 
where the first and second terms of Eq. ()17j) correspond to arrival at x = and x — 1, and 
where 

POO 

/ duuG(l,v,u) = 1 . (18) 

Integral equation (fTB^) follows immediately from this interpretation. The first, second, . . . 
terms in the iterative solution of the integral equation represent the mean time to reach the 
boundary for the first time, the mean time between the first and second boundary collisions, 
. . . , weighted with a factor p for each reflection. 

In Eq. (|16J) the asymptotic form of T(0, v) for small and large v is determined by the 
first and second terms, respectively, of the kernel (fTTjl . For p > the particle is not always 
absorbed the first time it reaches the boundary, which implies T(0,v) > Ti(0,v). For small 
v we look for a solution with the asymptotic form 

T(Q,v)~ v y , v\0, (19) 

with exponent j(p,r) smaller than the value 7 = \ for p = in Eq. (fTTj) . Substituting Eq. 
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(I19|) in Eq. ()16|) . we find that the asymptotic form is consistent with the integral equation if 

pr 7 = 2 sin ( 1 - 27 n ) . (20) 
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Solving Eq. (j2*Uj) for 7 graphically, one finds a unique solution satisfying < j(p, r) < \ 
in the physical region 0<r<l,0<p<l. Comparing Eqs. (JD) and (|2U|). one sees that 
7(p> r ) = 2</>(p, r). Thus, the exponents of v in the results (JIJ), (0) for the half line and (fTHj) . 
f!20|) for the finite line are the same. This suggests that the small v behavior (|19p. (|20|) is 
due to repeated collisions with the same boundary. 

For large v, T(0, v) ~ f ~ 1 + T(l, t>), corresponding to ballistic propagation from x = to 
a; = 1 in the time v . Combining this relation with Eqs. (jEJ) and ()12j) yields 

T(0,v) « + P T{0,rv) , u -»• 00 . (21) 

Equation (j21j) also follows from integral equation |T6|) on using the asymptotic forms fTTj) 
for large v and \F 2 (l; |, I; z) ~ iTr 1 / 2 ^ -1 / 4 exp(2z 1 / 2 ) for large z. For p < 1, the mean 
absorption time is expected to decrease to zero in the large v limit. Iteration of Eq. (|2T|) 
and/or substitution of the ansatz T(0,v) ~ _Bt>~^ in Eq. (|2*Tj) leads to the asymptotic forms 

{ r -1 
v , p < r < 1 , 
r- p * - ( 22 ) 
Bv -\np/\nr ^ r < p <1 , 

for v — > 00, with p < 1. 

The large w behavior for p = 1, r < 1, corresponding to inelastic reflection with zero 
absorption probability, requires special attention and is considered in Section V-C. 



III. NUMERICAL SOLUTION OF THE INTEGRAL EQUATION 

To interpolate between the asymptotic forms of T(0, v) for small and large v, given in the 
preceding Section, we solved integral equation numerically. Changing the integration 
variable from u to ru, we solved the resulting equation, 

T(Q,v) =Ti(0,u) + -^ / duuG{l,v,r- l u)T{<d,u) , v > , (23) 
^ Jo 

by iteration. 
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To evaluate the integral in Eq. ()23j) numerically, we changed the integration variable 
from u to 

z ~ g r 3 ' < u < B , 



w 



,3 



-Q: 



B < u < oo 



(24) 



g7s 

where a is a positive constant, and then used Simpson's rule, exact for polynomials of degree 
3, with equally spaced mesh points in z and w. To ensure good accuracy for small z (small 
u), we subtracted off the leading asymptotic form of the integrand, fit to a power law, and 
then integrated it analytically. In the region where both v and u are small, where the first 
term in the kernel ()17|) diverges, the integrand was fit to (v 3 + w 3 ) -1 times a power law in 
z and then integrated analytically. The parameters B and a and the mesh size were chosen 
to give good convergence with about 100 mesh points in each of the intervals < B < u 
and B < u < 00. Starting with Ti(0,t>) as the first approximation to T(0,t>), one typically 
needs about 10 iterations for convergence. 

IV. SIMULATIONS 

In our simulations the motion of the particle is governed by the difference equations 

/A 3 x \ 1/2 

+ V3r n+1 ) , (25) 
v n+ i = v n + (2A n+ i) 1/2 r n+ i . (26) 

Here x n and v n are the position and velocity at time t n , and A n+ i = t n+ i—t n . The quantities 
r n and s n are independent Gaussian random numbers such that 



(r n ) = (s n ) = 0, (rl) = (sl) = l. (27) 



As discussed in 



[16]. 



this algorithm generates trajectories which are consistent with the 
exact probability distribution P(x,v,t) of a randomly accelerated particle in free space, i.e. 
in the absence of boundaries. In free space there is no time step error in the algorithm. 
The time step A n+1 is arbitrary. However, close to the boundaries trajectories are not 
generated with the correct probability, because the free space distribution P(x, v, t) includes 
trajectories which wander outside the interval < x < 1 and return during the time t. 
As in we make the time step smaller near the boundaries to exclude these spurious 
trajectories. 
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Using the algorithm (|25j ) -([27 p with large time steps away from the boundaries, instead of 
a conventional algorithm with a constant time step, enables us to simulate the particle for 
much longer times. 

For a particle in free space with position and velocity x n , v n at time t n , the coordinate 
x n+ \ at a time A n+ i is distributed according to a Gaussian function jlfil ]. with a maximum 

1 /2 

at x n + v n A n+ i and root-mean square width (|A^ +1 ) . Choosing A n+ i so that 

x n + v n A n+1 - cA 3 n % > , (28) 
x n + v n A n+1 + cA% < 1 , (29) 

where the constant c is about 5 or larger, ensures that the Gaussian distribution lies well 
within the interval < x < 1. 

The largest A n+ i consistent with inequality (|2"£|) is given by D(x n ,v n ), where 

x + vD{x,v) -cD(x,v) 3/2 = . (30) 

To avoid solving Eq. ()30|) . cubic in D 1 ^ 2 , at each step of the algorithm, one may use the 
approximation 



D(x, v) 



x\ 2 / 3 ( v \ 2 
~c) + \c) ' ^>0 





> 2 V 




X 



(31) 
v < , 



which is asymptotically exact for v — > and v — > ±oo and accurate to better than 13% for 
— OO < V < oo. 

The most efficient time step to use in the algorithm is the largest A n+1 consistent with 
both inequalities (}2"%|) and (j2HJ). Using a smaller A n+ i slows the simulation without improving 
the accuracy. Since the inequality (|29|) follows from (J28)) on making the reflection symmetric 
substitution v n — > — v n , x n — > 1 — x n , the optimal time step is 

A n+ i = Mm{D(x n , v n ), D(l - x n , -v n )} + S , (32) 

where Min denotes the smaller of the two quantities inside the square brackets. As in 
a small minimum time step S has been included in Eq. (|32j) . Without it, the step size 
decreases to zero as the particle approaches the boundary, and, as in Zeno's paradox, the 
particle never gets there. After inclusion of 5, the particle not only reaches the boundary 
but jumps slightly pass it. We kept the overshoot small by choosing a sufficiently large value 
for the parameter c in Eqs. fl2"5j) - (f3"Tjh 
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In the case p = 1, r < r c of highly inelastic collisions with absorption probability zero, 
the speed of the particle becomes extremely small after many boundary collisions, and great 
care is required to simulate the behavior reliably According to Eqs. (|2T)j) and (j32)L the 
root-mean-square velocity change at each time step has the minimum average value 

Av = (25) 1/2 . (33) 

After each boundary collision we set Av equal to 1/500 of the velocity just after the collision. 
This Av is the smallest velocity the algorithm can reliably handle. The corresponding value 
of 5, given by Eq. is then used until the next boundary collision. 

The results for the mean time T(0, v) reported in Section V are based on 10 4 independent 
particle trajectories for each value of the initial velocity v. 



V. RESULTS 



A. Partially absorbing, elastic boundaries 

In this subsection we present our various results for p < 1, r = 1. In this regime the 
particle is absorbed at the boundary with probability 1 — p and reflected elastically with 
probability p. This is the "partial survival" model, studied on the half line in Refs. 

For p < 1, r = 1 the mean absorption time has the asymptotic forms ([19)1 . (|20|). and 
(|2*2*j) for small and large v , respectively. As the reflection probability p increases from to 
1, j(p) decreases from the Masoliver-Porra result 7(0) = ~ in Eq. (|TT|) for absorption on 
first arrival at the boundary, to the value 7(1) = 0, corresponding to no absorption at all, 
or T(0, v ) = 00. 

In Fig. ^our results for the mean absorption time T(0,v) for partially absorbing, elastic 
boundary conditions as a function of the initial velocity v are compared for r = 1 and 
p = 0.75, 0.5, 0.25, 0. The exact asymptotic forms (|19|) and (J2*U)| for small v and (|22|) for 
large v, the numerical solution of Eq. (|16J) . and the simulation results are clearly in excellent 
agreement. 

For fixed v the mean absorption time in Fig. ^ increases with increasing reflection prob- 
ability p, as expected. For fixed p the mean absorption time does not vary monotonically 
with v but has an absolute maximum at an intermediate velocity. For much larger initial 
velocities the particle is absorbed more rapidly since it bounces back and forth between 
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the boundaries, colliding at a rapid rate. For much smaller initial velocities the particle is 
absorbed more rapidly because it collides repeatedly with the boundary where it starts. 



B. Partially absorbing, inelastic boundaries 

Here we present our results for p < 1, r < 1. In this regime the particle is absorbed at the 
boundary with probability 1—p and reflected inelastically with probability p and coefficient 
of restitution r. Results for the half-line geometry with this boundary condition are given 
inRef. Q- 

In Fig. El our results for the mean absorption time T(0,t>) as a function of the initial 
velocity v are compared for (p,r) = (0.75,0.25), (0.501,0.5), and (0.25,0.75). Again the 
exact asymptotic forms ()19|) and (j20j) for small v and (}2*2*|) for large v , the numerical solution 
of Eq. (fT7)|) . and the simulation results are in excellent agreement. Note that the large v 
asymptotic forms ()22)1 for p < r < 1 and r < p < 1 are different and that the data test 
both. Corrections to the large v asymptotic form in the crossover regime p ~ r lead to 
conspicuously slower convergence for (p,r) = (0.501,0.5). 



C. Non-absorbing, inelastic boundaries 

In this subsection we consider boundaries with p — 1, r < 1. At the boundary the particle 
is reflected inelastically with probability 1 and absorbed with probability zero. This is the 
case of interest in connection with driven granular matter. Cornell et al. [5] predicted that 
for r < r c , where 

r c = e ^ /V ^ = 0.163..., (34) 

the particle undergoes "inelastic collapse," making an infinite number of boundary collisions 
in a finite time, coming to rest, and remaining there. 

The prediction that the inelastic collisions localize the particle at the boundary for r < r c 
was questioned by Florencio et al. 17. | and Anton on the basis of simulations. 

The equilibrium distribution function P(x, v) of a particle moving on the finite line be- 
tween two inelastic boundaries was studied by Burkhardt, Franklin, and Gawronski |lj| for 
r > r c and Burkhardt and Kotsev [l4| for r < r c with exact analytical and numerical cal- 
culations and simulations, similar to this paper. According to P(x,v) varies smoothly 
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and analytically with r throughout the interval < r < 1 and does not collapse onto the 
boundaries. For r < r c the equilibrium boundary collision rate is infinite, but the collisions 
do not localize the particle at the boundary (see footnote [loj]). 

In the case p = 1, r < 1 iterating integral equation (fTT)|) adds the mean time to reach the 
boundary for the first time, the mean time between the first and second boundary collisions, 
etc., taking the change in speed in each collision into account. Thus the solution T(0,v) 
represents the mean time T(0, v) for an infinite number of boundary collisions. 

For p = 1, r < r c , T(0,v) has the asymptotic forms ()19j) . (|20|) for small v. As the 
coefficient of restitution r increases from to 1, 7(7") decreases from the Masoliver-Porra 
result 7(0) = \ in Eq. (fTT|) to the value j(r c ) = 0, confirming the result of Cornell et al. for 
r c in Eq. The negative, unphysical value of 7 for r > r c signals the breakdown of the 

solution with finite T(0, v). For r > r c the mean time for an infinite number of collisions is 
infinite. 

In contrast with the results for p < 1 given in Eq. (}22|) . in the case p = 1, r < 1, the 
mean time T(0,v) does not vanish in the limit v — > 00. This may be understood by reverse 
iteration of the large v recurrence relation 1)21)1 . which yields 

T (0, r _n i>) « (r + r 2 + . . . + r") iT 1 + T(0, u) . (35) 

This result holds for any finite velocity v large enough so that the time to travel between 
the two boundaries is accurately given by the ballistic time v~ l . The first n terms on the 
right-hand side give the time a particle with initial velocity r~ n v takes to reach velocity v 
and sum to r(l — r)~ l v ~ l in the limit n — ► 00. 

According to Eq. 1)35)1 . lim^oo T (0, r _n f ) = T(0, cxd) is finite and non- vanishing for 
r < r c . The most general asymptotic form for large v consistent with these two properties 
and with the recurrence relation ()21jl is 

{00 , p = 1 r > r c , 

r , (36) 
/per(lnv) - w ' P = 1 , r < r c , 

where f peT (hiv) is a periodic function of Inv with period | lnr|. 

The periodic term in Eq. ()36)) came as a surprise to us. The large v recurrence relation 
()2*T|) implies T(0, f ) > T(0, rv) for p = 1, which suggests, but does not guarantee, that 
T(0, v) is a monotonically increasing function of v. The asymptotic form 1)36)1 for r < r c 
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with the periodic term is entirely consistent with T(0, v) > T(0, rv) and with Eq. (|35p. The 
periodicity already appears in the following crude approximation: For < v < v * define 
T(0,v) = Av 1 , in accordance with Eq. and then use the large t> recurrence relation, 

Eq. (J2U) with p = 1, to obtain T(0, u) for u > u*. 

To determine T(0, u) with simulations, we measured the mean time T/v(0, i>) for iV bound- 
ary collisions, plotted it versus iV -1 , and then extrapolated to iV — > oo. This is shown for 
r = 0.1 and v = 0.001 and 0.0001 in Fig. H 

In Fig. |3] our results for the mean time T(0, v) for an infinite number of boundary 
collisions are shown as a function of v for p = 1 and r = 0.1, 0.01, and 0.001. For fixed v 
the mean time increases with increasing r, as expected. 

Solving integral equation (jlfij) reliably becomes increasingly difficult as r decreases, and 
in Fig. El we only show the numerical solution for r = 0.1. The simulation data are in good 
agreement with this numerical solution, with the asymptotic form ()19j) . (|20|) for small v, and 
with the onset of periodic behavior in accordance with Eq. (JHHJ). The dashed lines for large 
v were calculated from the dashed lines for small v using recurrence relation (|2*T|) . 



VI. CONCLUDING REMARKS 



Masoliver and Porra calculated the mean absorption time T(x, v) of a randomly accel- 
erated particle on the finite line < x < 1, assuming that it is absorbed the first time it 
reaches either x = or x = 1. We have considered a more general boundary condition, 
parametrized by the reflection probability p and the coefficient of restitution r. We derived 
an integral equation which determines T(x,v) and used it to obtain the exact asymptotic 
form of T(0, v) for large and small v and numerical results for intermediate v. The asymp- 
totic forms and numerical results are in excellent agreement with our computer simulations 
of the randomly accelerated particle. 

The case p — 1, r < 1 corresponds to a single particle moving between walls with which it 
collides inelastically. This simple and fundamental system is of interest in connection with 
properties of driven granular media. We find that the mean time for an infinite number of 

)oundary collisions is finite for r < r c and infinite for r > r c , as predicted by Cornell et al. 

sj]. The prediction that the infinite sequence of collisions leads to inelastic collapse, with 
localization of the particle at the boundary, has been questioned on the basis of simulations 
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jl7, 18 1 and an analysis |lj] of the equilibrium distribution function P(x,v) using the same 
approach as in this paper. The results for the mean time reported here are largely inde- 
pendent of the inelastic collapse controversy. Both the integral equation and the simulation 
procedure sum the mean time, collision by collision, without directly addressing the question 
of localization at the boundary. 
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APPENDIX A: SOLUTION OF EQ. © 



With the substitution T(x,v) = T(x, —v), Eq. (0) takes the form 

dx dv 2 

Q 



f(x,v) = l. (Al) 



This differential equation is the same as Eq. (3) of Ref. j 1 31] - except for the inhomogeneous 
term 1 on the right-hand side. We follow Ref. [13j very closely in solving it. 
First we introduce the Laplace transform 



Q{s,v) = / dxe~ sx T{x,v) , (A2) 



o 

which according to Eq. (jAlj) satisfies 

d 2 \ ~ 1 
sv -q^)Q( s ' v ) = - + vT & v ) ■ ( A3 ) 

This the same as Eq. (A2) in except for the extra term s _1 on the right-hand side. As 
in Q], we solve Eq. (jA3|) in terms of Airy functions and invert the Laplace transform, with 
the help of the Faltung theorem. This yields 

1 r x / ,,3 



f(x,v) = ^M 1/2 / o dyy~ 1/2 e- v /18y I^ 1/6 lj 



,,1/2 roc /o,,3/2„,3/2\ 

V — jf duu^e-^^I^ {^—) T(0, u) 

e -v 3 /9(x-y) gfl q\ 

dy- oi { y,u) ^ v>Q ^ 



3V3r(|) Jo y (x - y) 2 /3 dv 
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which is the same as Eq. (6) of ^j], except for an additional term (the first term) on the 
right-hand side. 

To express the unknown function dT(y, 0) / dv on the right-hand side of Eq. (|A4J) in terms 
of the other unknown T(0,u), we first take the limit v = in Eq. (|A4j) . which yields 



f(x,0) 



3V3ro 



o2/3 poo 

— r(|)x 2 / 3 + X" 2 / 3 / duue- u:i / gx f(0,u) 
2 " Jo 

dy df(y,0) 



(A5) 



! (x — y) 2 / 3 dv 

From reflection symmetry T(x, 0) — T(l — a;, 0) = 0, and dT(y, 0)/dv = —dT(l — y, 0)/dv. 
Substituting Eq. ()A5|) in the first of these relations and making use of the second, we obtain 



dy df(y,0) 



/ \x — y| 2 / 3 dv 
where we have used 



Equation 
derived in 



duu 



e -u 3 /9x e -u 3 /9(l-x) 



x 2/3 (! _ x )2/3 



±w 2 + T(0,w) 



duu' 



„—u 3 /9x 



X 



2/3 



3 2/3 r( l )x 2/3 



(A6) 



(A7) 



\6|) may be regarded as an integral equation for dT(y,0)/dv. The solution, 



13( using the approach of 10. 
df(y,0) 



2CJ, is 



where the function R(y,u) is given in Eq. (fT5j). 
Substituting Eq. (|A"8l) in (lA"4l) yields 

T(x, f 1 



duw u) - i?(l - y, u)] \u 2 + T(0, it 



(A8) 



M = -^vf' 2 dyy-^e-^I^I^/A^ 



2-3V3r( 2 



e -t) 3 /9(a;-j/) roc 

~, ToTT / duu 3 [R(y, u) — R(l — y, u)] 



+ 



3; jo (x - y) 2/3 Jo 
duuG(x,v,u)T(0,u) , 



(A9) 



where G(x,v,u) is given in Eqs. ([T4"jl . Replacing T(x, u) by T(x,—v) in Eq. ()A9j) and 
evaluating the integral f °° duu 3 [. . .], we obtain the Green's function solution for T(x, —v) 
in Eqs. (fT3^) - (fT3j) and reproduce the Masoliver-Porra result for Tx(x,v) in Eq. (JTUJl . 



15 



[1] H. P. McKean, J. Math. Kyoto Univ. 2, 227 (1963). 
[2] M. Goldman, Ann. Math. Stat. 42 2150 (1971). 
[3] Y. G. Sinai, Theor. Math. Phys. 90, 219 (1992). 
[4] T. W. Burkhardt, J. Phys. A 26, L1157 (1993). 

[5] S. J. Cornell, M. R. Swift, A. J. Bray, Phys. Rev. Lett. 81, 1142 (1998). 
[6] M. R. Swift and A. J. Bray, Phys. Rev. E 59, R4721 (1999). 
[7] T. W. Burkhardt, J. Phys. A 33, L429 (2000); Phys. Rev E 63, 011111 (2001). 
[8] G. De Smedt, C. Godreche, and J. M. Luck, Europhys. Lett. 53, 438 (2001). 
[9] T. W. Burkhardt, J. Phys. A 30, L167 (1997). 
[10] J. Masoliver and J. M. Porra, Phys. Rev. Lett. 75, 189 (1995); Phys. Rev. E 53, 2243 (1996). 
[11] G. Fichera, in Boundary Problems in Differential Equations, edited by R. E. Langer (Univer- 
sity of Wisconsin Press, Madison, 1960); O. A. Oleinik, Sov. Math. 5, 1129 (1964); J. N. Yang 
and M. Shinozuka, J. Acoust. Soc. Am. 47, 393 (1970). 
[12] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products (Academic, New 
York, 1980). 

[13] T. W. Burkhardt, J. Franklin, and R. R. Gawronski, Phys. Rev. E 61, 2376 (2000). 
[14] T. W. Burkhardt and S. N. Kotsev, Phys. Rev. E 70, 026105 (2004). 

[15, E, m 8 tt» ^ GM _ 3(2„)-V/V/s („3 + ,,)- for a particle m „vi„ g on 

the half line x > with a single boundary at x = 0, derived by McKean [1J and independently 
by Cornell et al. |5J]. 

[16] D. J. Bicout and T. W. Burkhardt, J. Phys. A 33, 6835 (2000). 

[17] J. Florencio, F. C. Sa Barreto, and O. F. de Alcantara Bonfim, Phys. Rev. Lett. 84, 196 
(2000). 

[18] L. Anton, Phys. Rev. E 65, 047102 (2002). 

[19] This result, which may seem counter-intuitive, has the following mathematical origin: The 
equilibrium distribution P(x, v) obtained in lj, 1J] varies asymptotically as P(0, v) ~ v~@^ 
for u\0, where the exponent (3{r) increases monotonically from to | as r decreases from 
1 to 0, with (3{r c ) = 2. Thus, the equilibrium collision rate I = dvvP(0, v) diverges at 
the lower limit for r < r c . This result only holds for the continuous time dynamics (^Q). For 



16 



discrete dynamics with a nonvanishing time step, the number of collisions in a finite time is 
necessarily finite. 

[20] D. Porter and D. S. G. Sterling Integral Equations, (Cambridge University Press, Cambridge, 
1990). 



17 




FIG. 1: Mean absorption time versus initial velocity for r = 1 and, from top to bottom, p = 0.75, 
0.50, 0.25, and 0. The solid line for p = is the exact result of Masoliver and Porra. The other 
solid lines are the numerical solutions of the integral equation (|16|) . The points are the results of 
our simulations. The dashed lines show the exact asymptotic forms (|19|) . (|20[) for small v and (|22|) 
for large v. The proportionality constant in Eq. (|19|) was chosen to fit the simulation data. 
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FIG. 2: Mean absorption time versus initial velocity for, from top to bottom, (p, r) = (0.75,0.25), 
(0.501,0.5), (0.25,0.75), (0,r). The bottom curve is the exact result ()10|) of Masoliver and Porra. 
The other solid lines are the numerical solutions of the integral equation (|16|) . The points are the 
results of our simulations. The dashed lines show the exact asymptotic forms (|19[). (|20[) for small 
v and (|22|) for large v. The proportionality constant in Eq. (|19|) was chosen to fit the simulation 
data. Note the slow convergence to the asymptotic form for large v for (0.501,0.5) in the crossover 
region p ~ r. 
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FIG. 3: Mean time Tn(0,v) for N boundary collisions versus N 1 for r = 0.01 and, from top to 
bottom, v =0.001 and 0.0001. The mean time T(0, v) for an infinite number of collisions is estimated 
by extrapolating to the vertical axis. This yields T(0, 0.001) = 0.108 and T(0, 0.0001) = 0.039. 
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FIG. 4: Mean time for an infinite number of boundary collisions for p = 1 and, from top to bottom, 
r = 0.1, 0.01, 0.001. The upper curve is the numerical solution of the integral equation (|16|) for 
r = 0.1. The points are the results of our simulations, obtained by extrapolation, as in Fig. |3J 
The bottom curve (no simulations) is the exact Masoliver-Porra result 1)10(1 . The dashed lines on 
the left show the exact asymptotic form ()19l) . (|20l) for small v. The proportionality constant in Eq. 
(J19|) was chosen to fit the simulation data. For large v, approximately periodic behavior in Inv, as 
in Eq. (|36|) . is expected. The dashed lines for large v in Fig. |I]were calculated from the dashed 
lines for small v using recurrence relation 1)21(1 . 



